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'~~' Abstract: We investigate vortex solutions of a charged scalar field in Einstein- 

CN Maxwell theory in 3+1 dimensions with the addition of an axionic coupling to the 

f-f^ Maxwell field. We show that the inclusion of such a term, together with a suitable 

^ potential for the axion field, can induce an effective Chern-Simons term on the 2+1 

t::^- dimensional boundary. We obtain numerical solutions of the equations of motion 

^ and find Maxwell-Chern-Simons like magnetic vortex configurations, where the mag- 

netic field profile varies with the size of the effective Chern-Simons coupling. The 
axion field has a non-trivial profile inside the AdS bulk and on the 2+1 dimensional 
r*- boundary but does not condense at spatial infinity. 



1. Introduction 

The subject of holographic superconductivity (HS) has received much attention in 
the past 2 years since it was initially proposed in [1], [2] (see [3] for a recent review), 
[4] - [43] It has provided a fascinating and novel approach to the study of phase 
transitions in certain 2+1 dimensional systems by embedding them holographically 
in the background of an asymptotically AdS black hole. Phase transitions in the 
boundary theory appear in this setting as scalar hair that develops below a certain 
critical temperature when a charged scalar field is added to Einstein-Maxwell theory 
which supports either an AdS Schwarzchild or Reissner-Nordstrom type black hole 
solution. This bulk scalar field is known to be dual to certain operators in the 
boundary CFT and the development of scalar hair outside the black hole corresponds 
to these operators 'condensing', signifying a phase transition. 

Strictly speaking, in standard HS [1], [2], the dual CFT has no dynamical gauge 
field, so the phase transition is one that appears to describe the change in the super- 
fluid phase rather than a superconducting phase. ^ 

In [5] a different approach was taken by allowing for a 2;-independent term in 
the expansion of A^{r,z) = a^{r) + zJ^{r) + 0{z^). Here (r, x) correspond to 2d 
polar coordinates of the plane at fixed z, where z is the conformal coordinate of 
the 'intetrior' of the AdS^ background, with z = corresponding to the 2+1 di- 
mensional boundary at infinite distance from the AdS black hole. The expression 
-dra^{r) is then identified with a boundary magnetic field which one can think of as 
the result of weakly gauging the global f/(l) symmetry of the boundary CFT. The 
resulting solutions are rather non-trivial and depending on the various assumptions 
and boundary conditions imposed, were shown to describe either magnetic droplets 
or magnetic vortices living on the boundary (the authors of [9] obtained similar so- 
lutions although using a different approach where a magnetic field is added by hand 
on the boundary ai z = and its flux constrained to be quantised and related to the 
quantised superfiuid vorticity). 

Holographic supefiuid vortices have also been studied in [18], [19]. Here the 
approach is to take a^{r) = so there is no magnetic field. Solutions to the bulk 
Einstein-Maxwell theory were found which on the CFT side describe instead super- 
fluid vortex configurations. In this solution the superfiuid velocity is related to the 
current J^ appearing in the expansion of the A^ component of the bulk f/(l) gauge 
field near the CFT boundary aX z = 0. 

The appearance of such magnetic vortices (or superfiuid vortices) on the 2+1 
dimensional boundary at asymptotic infinity are triggered when the effective tem- 
perature of the system is above the critical value for which the scalar field condenses. 
They have quantised magnetic flux (or vorticity) and are thus topologically stable 

^A very interesting recent paper [36] discusses how the standard model of a HS can be modified 
to include a dynamical gauge field in the boundary CFT. 
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and are very reminiscent of the vortices discovered by Nielsen and Olesen in flat 2+1 
dimensional space-time [45]. 

There exist, however, different kinds of magnetic vortices in flat 2+1 dimensional 
space-time if one includes not just a Maxwell term but also a Chern-Simons (CS) 
term for the electromagnetic potential. Pure self-dual Chern-Simons vortices exist 
(when no Maxwell term is present) for suitably chosen charged scalar fleld potential 
[46], [47]. These pure CS vortex solutions are topologically stable, have quantised 
magnetic flux but in addition, unlike pure Maxwell type vortices, they also carry 
quantised electric flux. A distinguishing characteristic feature is a magnetic fleld that 
peaks in a ring outside of the 'core', unlike the pure Maxwell case where the magnetic 
field peaks inside the core. It is also possible to consider topological magnetic vortices 
in a theory which has both Maxwell and a Chern-Simons term present for the U{1) 
gauge field (see eg [48] for a comprehensive recent review). It is an interesting problem 
to see if such CS or mixed Maxwell-CS vortices can appear in a holographic setting. 

In the 'standard' model describing HS a complex scalar $ is coupled to a f/(l) 
gauge field described by a standard Maxwell action on an AdS4^ background. This 
scalar has a negative mass V{^) = —^\^\'^, where L is a scale defined by a non- 
vanishing cosmological constant A, which remains above the Breitenlohner-Freedman 
stability bound [49] and thus does not induce an instability of the system. 

The obvious problem one faces in trying to generalise the dynamics of the bulk 
gauge field to include a CS action is that no CS term exists in a 3+1 dimensional 
theory. ^ In principle one could add by hand a boundary action containing a pure 
CS term for the gauge field. However such an approach will lead to the appearance 
of explicit delta functions in the bulk equations of motion. Whilst the solution of this 
new system of equations may still be technically feasible, it is in some sense against 
the spirit of holography and we will not pursue this approach here. 

Instead we shall adopt a more indirect approach. As we will see, by including 
an additional pseudo-scalar axion field, 6, in our AdS^ background, with the usual 
coupling to the topological term F A F, we can induce an effective CS term for the 
f/(l) gauge fields that live on the 2+1 dimensional boundary. Such a term would 
then modify the magnetic vortex solutions found in [5], [9] and perhaps display 
properties of the more general fiat space Maxwell-CS vortices mentioned above such 
as the magnetic field peaking outside of the core. The coupling 6F A F can arise, for 
example, due to the axial anomaly if one included charge chiral fermions in the bulk 
theory. 

This paper is organised as follows: in section 2 we introduce the Maxwell the- 
ory with charged scalar in the AdS^ background, including, in addition, the axion 
field and its couplings. We make an ansatz for the dynamical fields, determine their 
equations of motion and discuss the relevance of terms in their near-boundary z = 

^gravity duals in 4+1 dimensions with both Maxweh and Chern-Simons terms present have been 
considered in [20], [21] 
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expansion. We shall see that the inclusion of the 6FAF term gives rise to an effective 
CS term on the boundary. Section 3 describes numerical solutions of this system, 
where we show that below a critical temperature the scalar field condenses and the 
corresponding magnetic field profile varies with the parameter k,, (the coefficient in 
front of the 6F A F term, which is related to the effective CS coupling on the bound- 
ary). In this section we also discuss the azimuthal current JAdSj fhe charge density 
PAds and the axion field 6{r) evaluated on the boundary and investigate the T > T^ 
region, where we show that the system admits a "trivial" solution where the scalar 
field vanishes asymptotically. In Section 4 we perform a near-boundary expansion of 
the equations of motion to try and uncover the asymptotic analytic properties of the 
numerical solutions for large and small distances in the 2+1 dimensional boundary . 
Finally, in Section 5, we provide a short discussion of the results and directions for 
further research. 



2. The Model 

The action we consider is: 

where 

A = -^,,c = -/.uF^" + ^,\D,^' - -^|$r + -^ee^'^'-^F.^F,^ + {d,ef + v{e) 

(2.2) 
where Gtv is the gravitational coupling, g the gauge coupling, D^ = d^ — iA^ 
and V{6) is the axion potential. The constant k, is, as we shall soon see, related to 
an effective CS coupling on the boundary. In our coordinate notation, z denotes a 
"radial" or "bulk" direction of the AdS4 space with z = being its boundary at 
infinity and (r, x) acting as coordinates for two-dimensional planes of constant z. 
The above form of the action in which the gauge coupling 1/g'^ appears only as an 
overall multiplicative factor in front of the 'matter' Lagrangian is useful in working 
in the probe limit (which we shall assume) where g is taken very large. To bring the 
action into this form, the canonical fields in the Lagrangian are redefined as 

A^^-A^, <^^-<^, e^-9, K^gn (2.3) 

9 9 9 

These rescalings allow us to work in the formal decoupling limit of (7 — ?■ oo 
where one can safely assume that the background space is fixed, uncharged and 
suffers from no back-reaction from fields in the Lagrangian. The space time is the 
A(iS'4-Schwarzschild black hole described by the metric 
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ds' = ^{-f{z)df + dr' + r'dx') + Sr-,dz^ 

Z^ Z-^j'yZ) 

f(z) = 1 - z' 

v^ = ^ (2.4) 

with z = 1 the location of the black hole horizon. This has a Gibbons-Hawking 
[44] temperature T, which is dual to the CFT temperature T via the AdS/CFT 
dictionary. For the purpose of a HS one is interested in the dependence of a dimen- 
sionless measure of the condensate as a function of a dimensionless measure of the 
temperature (as explaine din [3] ), which we shall discuss later in Section 3. 

We work with an ansatz of the form 
A, = 

Ar = 

$ = 0(r, 2)e™^ 
A^ = A^{r,z) 
Ao = Ao{r,z) 

e = e{r,z) 

V{e) = mp^ + \e, \/($) = --^|<l>|2 (2.5) 

Note in the above we have included the possibility of a linear term on the axion 
potential. The presence of the axionic coupling 6F A F already violates 6 — )■ —6 
parity so in principle such a term could be included. However as we shall see later, 
the coefficient A of the linear term must be related to the value of 6 evaluated on the 
boundary at 2; = in order to solve the equations of motion. This necessarily means 
that taking X ^ 0, 6 must be a constant everywhere on the boundary and so has a 
trivial profile there. It is important to point out that this constant value of 6 should 
not be thought of as the axion field condensing in this system. This is because any 
constant appearing in 6 can be effectively removed by adding an explicit term F AF 
to the bulk action with an appropriate coefficient. Thus no axionic "hair" exists. 
This is consistent with general arguments (see e.g. [3]) that the ^^5*4 black hole can 
only develop hair from charged scalar fields. ^ 



^Note that black hole solutions in 3+1 dimensions are known to exist which do indeed develop 
'conventional' axionic hair (by conventional we are excluding the possibility of stringy axionic hair 
that can arise when one has a dilaton and Kalb-Ramond field present which arise from a string 
theory compactification) when one includes the term 9F A F [50], [51]. However such hair only 
arises when one includes torsion into the background geometry. 
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With the above ansatz we obtain the following bulk equations of motion 



^4^2 I d^ifdzA^) + rdr [ -drA^ ] - -^<P\A^ -n)]- SKe^'^^^F.^d^e = 



z'd, ^d,6 + -dJrdrd)) + 



^'d. (^9.^) + h{rm - z^>.FF +~ = 0, (2.6) 

In what follows we will be interested in the behaviour of the various fields near 
the boundary z = 0. We shall assume the fields have the following expansions for 
small z 



(t){r,z) = z(l)o{r) + z^cpi^r) 
Ao{r, z) = fi + zp + z^A2 + z'-^As + z^A^. + .. 
A^{r, z) = a^ + zJ^ + z'^D2 + .... 

9{r, z) = 9o{r) + zdi{r) + z^^aW + ... (2.7) 

The AdS/CFT correspondence [52] teaches us that terms in the near boundary 
expansion of bulk fields act as sources for conformal operators in the dual theory. 
More precisely, 0i sets the vacuum expectation value of an operator Oi with dimen- 
sion i. For consistency, only one of these operators can exist at one time so focussing 
on the 01 case will exclude 02 and vice-verse. We will only consider 0i in the follow- 
ing. Within the context of HS the scalar field describes, in the dual CFT, the order 
parameter for the spontaneous breaking of the U{1) symmetry. Similarly the first 
and second terms in the expansion of Aq are associated to the Maxwell contribution 
to the chemical potential and charge density respectively, whilst those of A-^ are re- 
lated to the magnetic field B = -dra^ and azimuthal current density J^. With the 
presence of a Chern-Simons boundary term the precise definition of these quantities 
acquires a further contribution. For the detailed relations we refer the reader to the 
Appendix. 

3. Solution 

We are interested in solving the coupled partial differential equations eq(2.6). In 
order to avoid divergences at r = we make the following field redefinitions 

= zr"'0 
A = r'A^, (3.1) 
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such that, near the boundary, 



= 00 + #1 (3.2) 

A^ = Aq + zAi + z^A2 + ... (3.3) 

This imphes that the magnetic field becomes 

B = 2Aq + rdrAo. (3.4) 

and the condensate is 

0i(^) = 2^'0 = ^>i- (3.5) 

With these redefinitions and setting L = 1, the equations of motion become 






X ^2 



:-9,^9,Ao + 9,^9,Ao) = (3.7) 



rf 



d^Ao + - ( d^.Ao + -drA, ] - —02^0 



1 / 1 \ r^" 

- (^9,^0 + -drAoj - J 

+ (-r^S.^S.A^ + r^d,9drA^ + 2r9,^A;,) = (3.^ 






Z^ 



iKZ 



2 



+ ^^(-r2(9^Ao(9^4 + ^^-^^^-^^^o + 2rA^d-,Ao) = (3.9) 

We are interested in solving these equations to extract the boundary behaviour 
of the fields. Unsurprisingly the equations cannot be solved analytically, they do 
however present numerical solutions. 

Before investigating these, an important observation to make is that by compar- 
ing the Aq or A^ equations of motion with those derived from a purely 2+1 dimen- 
sional Maxwell-CS theory with charged scalar matter, we see the term ndzO ~ K,6i{r) 
(whose origin is the axion coupling 6F A F ) generates an effective CS coupling in 
the 2+1 dimensional boundary theory at 2; = 0. 
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3.1 Solution 1 

For the numerical procedure, we must specify boundary conditions on the fields. Ini- 
tially, we work with /i = const, J^ = const, 60 = I and n = 1. Then, for consistency 
at the boundary z = 0, one finds m^ = 1 and A = —2. In all solutions found, the 
scalar field (Figure 1) condenses independently of the parameter k. This is in close 
agreement with the results of [9] in which the scalar field vanishes at the location of 
the vortex. In Figure 2 we present the results for the magnetic field B{r). We find 
that for low enough k the system is dominated by a pure Maxwell solution, where 
the magnetic field peaks at the origin and vanishes for large r [5]. However, as k 
is increased the solution is driven away from the pure Maxwell configuration until 
eventually, at a critical value of k, the magnetic field vanishes at the origin and peaks 
at finite distance from it (red line), once again decaying to zero at large r. This 
behaviour is characteristic of self-dual Chern-Simons vortices in flat space [47] where 
B{r) forms a ring surrounding the origin. At even higher values of k, the magnetic 
field becomes negative close to the origin. The reader might at this point be worried 
about flux quantization, however the flux $b = J rBdrdx always remains positive 
due to uneven contributions to the integral from below and above the zero axis, hence 
we expect to maintain flux quantisation as required for topological solitonic solutions. 
The Charge density, which depends on the parameter k, is plotted in Figure 3. It 
tends to a constant value at asymptotically large values of r approximately indepen- 
dent of the value of k, at least for the values chosen. The effect closer to the core of 
the vortex, which appears to be a depletion of electric charge is in fact a temperature 
dependent phenomenon. This was flrst found in [5] where the authors showed that 
upon lowering the temperature in the region T < T^ the charge density switches 
from increasing to decreasing there. This effect, if not simply due to the complex 
numerical procedure, was suggested to be indicative of strong screening effects from 
the core of the vortex. However one should not interpret this as signifying that the 
vortex cannot support electric charge, as indeed one expects a Chern-Simons vortex 
to do. 

It is also interesting to study the behaviour of the axion fleld 6i{r). In Figure 4 
we show plots for different values of k,. Increasing the parameter k, pushes the axionic 
field to lower values whilst preserving its shape. The importance of 6i{r) is that for 
large r it plays the role of a CS coupling on the boundary, as discussed earlier. As 
we see from Figure 4, 61 tends to a constant for large r, which lends further weight 
to this interpretation. 

In passing we also mention there exists another set of solutions for which, rather 
than imposing Dirichlet conditions on fi, we allow for this to be dynamically de- 
termined by the numerical procedure whilst keeping p = const. Then by having 
a non-constant p at the boundary one can observe a pure Chern-Simons contribu- 
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tion to the azimuthal current J^, which was previously determined only by it's pure 
Maxwell component set manually to a constant (see Appendix). We leave the nu- 
merical exploration of such solutions, and the corresponding profile for the modified 
azimuthal current to future work [53]. 

3.2 Solution 2 

This section is devoted to analysing the set of solutions for T > T^, the critical 
temperature at which the scalar field condenses. As explained in [3] we are interested 
in the dependence of the condensate as a function of a dimensionless measure of the 
temperature. We take - for this measure, where T is simply the Gibbons-Hawking 
Temperature of the black hole and remains constant for varying parameters of the 
numerical procedure. Hence, varying values of /i allow us to change the effective 
temperature of the CFT. Figure 5 shows the profile of the condensate for different 
values of fi. As we make /i smaller, and hence the effective Temperature larger, 
the solution enters the region oi T > T^ and the scalar field doesn't condense. For 
increasingly small fi, the solution is driven towards 0i = 0. Conversely, when one 
raises fi, the profile of the solution remains that of the condensing scalar (blue line) 
as the system is well within the T < Tc region. 

4. Asymptotic analysis 

Even though the set of coupled equations eq(3.6) is highly non- linear and can only 
be solved numerically, one can try and perform an asymptotic analysis to get an 
analytical feel for the behaviour of the functions at the z = boundary. In the 
case of purely 2+1 dimensional CS vortices or mixed Maxwell-CS vortices, such 
an analysis proves useful in providing analytic expressions for the various fields in 
the regions near the core and far away from the vortex (see e.g. [48] for a detailed 
analysis). These analytic expressions, apart from their own intrinsic interest, can 
provide further checks on any full numeric solution. 

In our model we first need to obtain a series of equations which are derived from 
the full equations of motion by taking a small z expansion. For the origin of the 
various fields appearing in this expansion, we refer the reader to eq(2.7). From the 
equation of motion of the scalar field, expanding the fields around z = using eq(3.2) 
and eq(2.7), at 0{z) we find 

2nAo4>i-r^Al4>i + (l + 2n)^^ = 0. (4.1) 

r 

Similarly, the Aq equation gives 

drp(^ + IGrKeA + d^.p + lQKrpe2 + A3 = 0, (4.2) 



at order 0{z), and 

+16fi:rii(4^2 + rS^^a) + rA^ + 16/tr^3(2io + rAi) = (4.3) 

at order 0{z^). These equations govern the behaviour of the magnetic field, the 
condensate and the charge density shown in the solutions above. They are a com- 
plex system of coupled partial differential equations. In particular, with a simple 
rearrangement eq(4.1) can be expressed as 




^o = -2 + ^2\jy + (^ + ^^)^^^ (4-4) 

which implies that the magnetic field eq(3.4) 

^ ^ il + 2n){-r{dr4>iy + Mdr{rdr^i))) ^^ ^^ 

Then, to observe the full dependence of the magnetic field on k, one would have 
to substitute eq(4.3) rearranged for 0i in eq(4.5), even for k 7^ this is easy to do 
but, for the sake of simplicity, we will not quote the full result. Instead, for the 
simplifying choice of k = 0, where one has pure Maxwell, note that one can read off 
the asymptotic form of the condensate, and eq(4.3) becomes 



, , . ni dr{rdrD2{r))+rA,{r) 
0i(r)=r0i = W — , (4.6) 

so the condensate contains functions appearing in the expansion of A^ and Aq 
(recall that in our numerical simulation /x is kept constant). 

Indeed the lesson we are learning is that the set of equations that arise at each 
order in z do not form a closed set, in that the number of functions is always greater 
than the number of equations. Even at lowest order we see from the above that addi- 
tional 'source' terms D2{r), A4^{r) appear in the 3 equations that determine A{r), B{r) 
and 01 (r) and from which one could try and extract the small and large r behaviour. 
But to do so requires some knowledge of the behaviour of the functions D2{r), A^lr) 
in these same limits. 

Thus the situation is, unfortunately rather more complex than in the case of 
strictly 2+1 dimensional Maxwell-CS vortices [48]. One way of proceeding is to 
make approximate analytic fits to the functions D2{r) and A^lr) from their numerical 
profiles (which can be derived though we have not presented the plots here) for small 
and large r. Substituting these into eq(4.6) would then allow one to extract analytic 
asymptotic expressions for the condensate (pi, B{r) and A{r). In doing this one 
would impose the appropriate boundary conditions on the latter fields for r — )■ and 
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r — )■ oo in order to to yield finite energy solutions on the 2+1 dimensional boundary 
ai z = 0. We will leave this analysis to future work [53]. 

5. Discussion 

In this paper we have made a first attempt to find a holographic description of 
magnetic vortices in 2+1 dimensions where the the abelian gauge field is governed 
by both a Maxwell term and an effective CS term. We have achieved this within a 3+1 
dimensional theory describing an AdS black hole (in the probe limit), by including 
an additional neutral scalar field, the axion, and its coupling to the topological term 
F A F. The resulting equations of motion are, like the original HS, rather difficult 
to solve analytically and we found instead numerical solutions for the magnetic field, 
scalar condensate, charge density and axion fields as functions of the coordinates 
{r,z). The corresponding profiles of these fields on the boundary at 2; = were 
presented. For a range of values of the chemical potential fi corresponding to T < T^ 
they clearly describe magnetic vortex like solutions. An important difference when 
compared to those found in the literature [5], [9], is that the magnetic vortices we 
have found have the magnetic field peaking outside of the core when the parameter 
K, that multiplies the 6F A F term in the action, (and which is related to the effective 
CS coupling on the boundary), is above a certain value. Such a property of the 
magnetic field 'peaking' in a ring outside the core is very reminiscent of the profile 
of a magnetic field outside of a CS vortex [46], [47]. 

By lowering the values of fi we have also explored the region T > T^ and have 
found numerical profiles that show the condensate gradually vanishing at large dis- 
tances as one decreases /x. 

A further novel behaviour we have observed is the field reversal of the magnetic 
field near to the vortex core (in the case T < T^) as one continues to increase the 
parameter k. Whilst the phenomenon of distance dependent field reversal in magnetic 
vortices is known to occur in anisotropic superconductors [54] , its appearance in our 
model and its dependence on k would benefit from further investigation. One clue 
could come from the observation made earlier, that it is k9i that plays the role of an 
effective CS coupling on the 2+1 dimensional boundary. But 9i is a function of r, is 
dynamically determined and implicitly depends on the value of /t. It only tends to a 
constant at large r as can be seen from Figure 4. So in comparing with the analysis 
of standard 2+1 dimensional Maxwell-CS vortices we have to bear in mind that we 
have, effectively, an r dependent CS coupling that is only constant for large r. 

Another context in which oscillatory behaviour of a magnetic field is observed 
is when considering non-relativistic vortices in external magnetic and electric fields. 
Such theories are relevant in the ZHK model of the quantum Hall effect [55] where 
vortices occur in a statistical gauge field governed by a pure CS action. The magnetic 
CS vortices in this context can exhibit oscillatory or 'field reversal' behaviour similar 
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to what we have seen in Figure 2 (see section 4 of [48] for a review). This system 
is of course different from the one we have considered in this paper, since there is 
a distinction between the fictitious or 'statistical' CS gauge field and the physical 
background magnetic field. It would be worthwhile to investigate how such non- 
relativistic CS vortices could arise in a holographic framework. 
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7. Appendix 

The AdS/CFT dictionary states that each term in the expansion of the fields at 
2; — )■ is a source for a conformal field at the boundary. More precisely, for a 
conformal field < Oi> sourced by terms in the expansion of Ai 

where x are boundary coordinates in the dual theory, Son-sheii is the on shell, 
Euclidean version of the action eq(2.2), /3 is the inverse temperature and V is the 
spatial volume of the dual theory. The calculation to derive the Maxwell contribution 
to the conformal fields is illustrated in [5] . In our case we have in addition the 6F A F 
term which contributes a CS like term 29oez^paA^Fp^ when one considers boundary 
contributions to the variation 5 Son-sheii- 

Hence, including variations of this term with respect to progressive orders in the 
expansion of Ai at the boundary 2 = determine the full conformal fields. This 
gives. 



pAds = -o^^^o + '^OodrA^ 



Jlds = 2^-^x - 4^o9.Ao 



(7.2) 



where importantly we see the non-vanishing contribution of the effective CS 
terms in the definition of the charge density and azimuthal current. In the above 
we have ignored the contribution for the black hole horizon at 2; = 1 in evaluating 

OOon— shell- 
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Figure 1: The Condensate 4>i{r) = r(/>i(r) (in units of ;U = 1) for different values of the 
parameter k. As expected, the condensate vanishes where the vortex is placed. There is 
negligible variation in the condensate profile for this range of k. 
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Magnetic Field B(r) 




Figure 2: The Magnetic Field B{r) for different values of the parameter k. For low values 
of K with the magnetic field peaks at the origin (purple line) as per the pure Maxwell case. 
As one increases n this profile is modified (green line) until at a critical n = 0.6055 we see 
that the Magnetic Field vanishes at the origin and forms a ring around it (red line). For 
increasing k the magnetic field falls below the zero axis. 
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The Charge Densit?/ p(^)r^ 




Figure 3: The Charge Density pAdsi^) for different values of the parameter k. The charge 
density tends to a constant for large r and shows little modification with varying values of 

K. 
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Figure 4: The Axion Field Oi{r) for different values of k. The Axion Field is negative 
and increases with smaller values of k. 
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Figure 5: The condensate (pi^r) (in units of /x = 1) for different values of the chemical 
potential fi, and hence different temperatures. One sees that below Tc the scalar field 
condenses, but as one raises the temperature the scalar field drops gradually to zero. 
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